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Abstract. For a class of polynomial non-autonomous differential equations of 
degree n, we use phase plane analysis to show that each equation in this class 
has n periodic solutions. The result implies that certain rigid two-dimensional 
systems have at most one limit cycle which appears through multiple Hopf 
bifurcation. 



1. Introduction 
We consider differential equations of the form 

z:= — = z n + Px^z"- 1 + • • • + P n -x{t)z + P n (t) (1.1) 
at 

where z is a complex-valued function and P, are real-valued continuous functions. 
We denote by z(t, c) the solution of (|l.lfl satisfying z(Q, c) — c. Take a fixed real 
number u, we define the set Q to be the set of all complex numbers c such that 
z(t,c) is defined for all t in the interval [0,w]; the set Q is an open set. On Q we 
define the displacement function q by 

q(c) = z(u>, c) — c. 

Zeros of q identify initial points of solutions of which satisfy the boundary 

conditions z(0) = z(u>). We describe such solutions as periodic even when the 
functions Pi are not themselves periodic. However, if Pj are w-periodic then these 
solutions are also w-periodic. The main concern is to estimate the number of 
periodic solutions. This problem was suggested by C. Pugh as a version of Hilbert's 
sixteenth problem; it is listed by S. Smale as Problem 7 in [Tl|. Equations 
have been studied in detail by Lloyd in 1 1 1 j , using the methods of complex analysis 
and topological dynamics. 

Note that q is holomorphic on Q. The multiplicity of a periodic solution ip is that 
of <p(0) as a zero of q. It is useful to work with a complex dependent variable. The 
reason is that the number of zeros of a holomorphic function in a bounded region 
of the complex plane cannot be changed by small perturbations of the function. 
Hence, periodic solutions cannot then be destroyed by small perturbations of the 
right-hand side of the equation; periodic solutions can be created or destroyed only 
at infinity. Suppose that (p is a periodic solution of multiplicity k. By applying 
Rouche's theorem to the function q, for any sufficiently small perturbations of the 
equation, there are precisely k periodic solutions in a neighborhood of (p (counting 
multiplicity). Upper bounds on the number of periodic solutions of (|1 . 1|) can be 
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used as upper bounds on the number of periodic solutions when z is limited to be 
real-valued. This is the reason that Pi are not allowed to be complex- valued. 

When n = 3, equation is known as the Abel differential equation. It was 
shown in ^1] and J3j that Abel differential equation has exactly three periodic 
solutions provided account is taken of multiplicity. We describe equation Ijl.lfl as of 
Abel form. For n > 4, Lins Neto JOj has given examples which demonstrate that 
there is no upper bound, in terms of n only, on the number of periodic solutions. 
However, there are upper bounds for certain classes of equations. It was shown in 
that each of the following equations has n periodic solutions. 

z = z n + a(t)z, 

z = z n + a(t)z ]1 i 1 + {3(t)z, 

where n is odd in the second equation. Ilyashenko, in gave the following upper 
bound for the number of periodic solutions 

8 exp ((3C + 2) exp (|(2C + 3)")) 

where, C > 1 is an upper bound for the absolute values of the coefficients Pj(t). 
Although this estimate is non-realistic, it is the only known explicit estimate. In 
|12|. Panov considered the equation 

i = z n + a(t)z 2 + j3(t)z + j(t) 

and proved that the equation has at most three real periodic solutions if n is odd. 
Quartic equations having at least ten periodic solutions were described in p^; the 
coefficients were polynomial functions in t of degree 3. 
In this paper, we consider the class of equations 

z = z n + a(t)z n - 1 +P(t)z n - 2 (1.2) 

as a generalized Abel differential equation. We show that if (3(t) < 0, then this 
equation has at most two non-zero periodic solutions. We give conditions on a and 
[3 that imply the equation has exactly two non-zero periodic solutions, one non-zero 
periodic solution, or no non-zero periodic solutions. Particular cases of this result, 
with n = 4 and n = 5, were given in 0] and In Section 2, we describe the phase 
portrait of and recall some results from lT|. In Section 3, we state and prove 
our main result. In the final Section, we use the result of Section 3 to show that 
a certain family of rigid two-dimensional systems has at most one limit cycle and 
this limit cycle appears through multiple Hopf bifurcation. 

2. The Phase Portrait 

We identify equation with the n— tuples (Pi, P2, Pn) and write C for the 
set of all equations of this form. With the usual definitions of additions and scalar 
multiplications, C is a linear space; it is a normed space if for P = (Pi, Pi, ...,P n ) 
we define 

||P|| = max{ max |Pi(f)|, max |Pa(t)|,-> - , max |P n (t)|} 

0<t<uj 0<t<u Q<t<u 

The displacement function q is holomorphic on the open set Q. Moreover, q de- 
pends continuously on P with the above norm on C and the topology of uniform 
convergence on compact sets on the set of holomorphic functions. 
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The positive real axis and the negative real axis are invariant. Moreover, if ip is 
a non-real solution which is periodic, then so is ip, its complex conjugate. 

In JT], it was shown that the phase portrait of is as shown in Figure 1. 

We refer to for the details. There, the coefficients Pi(t) were 10— periodic. It 
can be verified that the same methods are applicable to the study of the number 
of solutions that satisfy z(0) = z(u>) whether the coefficients are periodic or not. 




Figure 1 . Phase Portrait 

Note that the radius, p, of the disc D is a sufficiently large number that depends 
only on ||P|| and u. If z = re 10 then the sets Gk, k = 0, 1, . . . , In — 3, which are the 
"arms" in the figure, are defined by 

_ r i kir a k-K a, 

G k = {z\r > p, - - < 6 < + -} 

n — 1 r rt — 1 r 

where a = max{6, 6||P||}. Between the arms are the sets Hk, k = 0, 1, . . . , 2n — 3, 
which are defined by 

kir a n (k + l)ir o, 

H k = {zr>p, +-<0<S 1 } 

n — 1 r n — 1 r 

For even k, trajectories can enter Gk only across r = p, and for odd k, trajectories 

can leave Gk only across r = p. No solution can become infinite in Hk as time either 

increases or decreases. Every solution enters D. Solutions become unbounded if 

and only if they remain in one of the arms Gk , tending to infinity as t increases if 

k is even and as t decreases if k is odd. 

Let q(P,c) = zp(uj,c) — c, where zp(t, c) is the solution of P 6 C satisfying 

zp(0,c) = c. Suppose that (Pj) and (cj) are sequences in L and C, respectively, 

such that q(Pj, Cj) — 0. If Pj P and cj — > c as j — > oo, then either g(P, c) = 0, 

in this case Zp(t, c) is a periodic solution, or zp(t,c) is not defined for the whole 

interval < t < to. In the later case, we say that zp(t,c) is a singular periodic 

solution. We also say that P has a singular periodic solution if Cj — > oo; in this case 

there are r and c such that the solution zp with zp(r) = c becomes unbounded in 

finite time as t increases and as t decreases. We summarize the results of 

Proposition 2.1. (i) Let A be the subset of C consisting of all equations which 
have no singular periodic solutions. The set A is open in C. All equations 
in the same components of A have the same number of periodic solutions. 
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(ii) The equation z = z n has exactly n periodic solutions. 

(iii) For each k and for each to and r > p, there is a unique 9o such that the 
solution z(t,c) with z(to,c) = c remains in G& for t > to (k is even) or 
t < to (k is odd) if and only if c = re 1 . 

3. Main Result 

Assume that P n (t) = 0. We call the solution z = 0a center if z(t, c) is periodic 
for all c in a neighborhood of 0. If the term z 11 in is replaced by Po{t)z n , 

then there are equations with a center when Pq has zeros. For cubic equations, this 
is related to the classical center problem of polynomial two-dimensional systems; 
we refer to |S] for details. However, when Po has no zeros then z — is never a 
center. Particular cases of this result were given in and 0] for n — 4 and n = 5, 
respectively. We give a brief proof, for the sake of completeness. 

Theorem 3.1. The solution z — is isolated as a periodic solution of with 
P„(i)=0. 

Proof. Suppose, if possible, that there is a open set U C C containing the origin 
such that all solutions starting in U are periodic. Then q = in the component of 
its domain of definition containing the origin. But the real zeros of q are contained 
in the disc D. Thus 

inf{c £ R : c > 0, z(t, c) is not defined for < t < u} < oo 

It follows that there is a real singular periodic solution; but a positive real periodic 
solution which tends to infinity can do so only as t increases. This is a contradiction, 
and the result follows. □ 

Now, we give the result about the number of periodic solutions. 

Theorem 3.2. Suppose that (3{t) < 0. Equation (|1.2(l has exactly n periodic solu- 
tions. 

Proof. With z = re l9 , we have 

6 = r n - x sin (n - 1)6 + r n ~ 2 a{t) sin (n - 2)6 + r n - 3 /3(t) sin (n - 3)0. 

If |c| > p and is real then the real solution z(t, c) remains outside the disk D either 
when t increases or when t decreases, and will become infinite. Solutions that enter 
Go or G n -i will leave Go or G„_i, except the solution that enters at pe l6 ° described 
in part (iii) of Proposition 12.11 this solution is real because any solution which is 
once real is always real. Therefore, the unique solution that becomes infinite is a 
real solution iffc = 0or/c = n— 1. On the other hand, no real solution is unbounded 
as t increases and decreases. Hence, no singular periodic solution enters Go or G„_i 
because singular periodic solutions are unbounded both as t increases and decreases. 
Thus, a singular periodic solution enters D from a Gk with odd k =/= n — l and leaves 
D to a Gj with even j ^ 0, n — 1. 

Now, if k is odd and 1 < k < n — 3, let 6\ = -^j^, 6 2 — I-^ ' anc ^ cons ider 

0(01 ) = r™- 1 sin ((n - l)0i) + r n - 2 a(t) sin ((n - 2)6 X ) + r n ~ 3 (3(t) sin ((n - 3)0i), 
0(0 2 ) = r n ~ l sin ((n - 1)0 2 ) + r n - 2 a{t) sin ((n - 2)6 2 ) + r n - 3 [3(t) sin ((n - 3)0 2 ) 
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But, 

sin ((n - 1)6*1) = sin (kn H ) < 0, 

n — 2 

sin ((n — 3)6*1) = sin (kn ) > 0, 

n — 2 

sin ((n - l)0 a ) = sin ((Jfe - 1W + ^— ^ ) > 0, 

n — 2 

sin ((n - 3)02) = sin ((A - 1W - ^ ~ ^f ) < 0, 

n — 2 

sin ((n - 2)6>i) = sin ((n - 2)6 2 ) = 0. 
Under the above hypotheses, 6(61) < and 8(62) > 0. Since p is a sufficiently large 
number, we assume that p > a (" 2 K n — 11 xhis condition on p guarantees that the 
arc of intersection of Gk with D lies inside the sector 

{r< P ,e 2 <e<e 1 }. 

Precisely, the condition implies that 

kn kit a (k — lW kit a 
> + -,- i— < 



n ~ 2 n — 1 p n — 2 n — 1 p 

Hence, solutions do not leave the sector (see Figure 2). Therefore, no singular 
periodic solution can enter D from Gk and leave D to Gk-i or Gfc+i. Since the 
phase portrait is symmetric about the x— axis, if k is odd and n < k < 2n — 3, no 
singular periodic solution can enter D from Gk and leaves D to Gfc_i or Gk+i- It 
follows that the equation does not have a singular periodic solution. 
Now, consider the class of equations 



sa 



with < s < 1. Since, s/3(t) < 0, any equation in this family does not have singular 
periodic solutions. The equation z = z n belongs to this family and has n periodic 
solutions. By part (i) of Proposition 12.11 each of these equations has n periodic 
solutions. □ 



Gk+i 




FIGURE 2. Sector around Gk, odd k 
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Lemma 3.3. Let A = f£ a(t)dt and B = f£ f3(t)dt. The solution z = of 
has multiplicity 

(i) n - 2 if B j= 

(ii) n — 1 if B = and A 0. 

(iii) n if A = B = 0. 

Proof. We write 

oo 

z(t,c) = ^a k {t)c k 

and substitute directly into the equation l|1.2|l . This gives a recursive set of linear 
differential equations for the a k (t) with initial conditions ai(0) = 1 and cife(O) = 
if fc > 1. The multiplicity is K if ai(u>) — 1, aj.(cj) = for 2 < k < K — 1 and 
ok(w) 0. Direct computations give 

ai(i) = l;ofc(t) = 0,2<fc<n-3, 

a«-2 = On-i = a;a n = 1. 

Solving the last three equations gives a n _2(w) = S, a n _i(a;) = A, and a n (u>) = uj. 
The result follows. □ 

Corollary 3.4. Consider equation (|1.2(l with f3(t) < 0. The equation has 

(i) two non-zero periodic solutions if B < 0; at most one is positive and at most 
one is negative. If B is small and negative, then there are two non-zero real 
periodic solutions, 

(ii) one non-zero periodic solution if B = and A ^ 0; this solution is a real 
solution. It is positive if A < and it is negative if A > 0, 

(iii) no non-zero periodic solutions if A = B = 0. 

Proof. Consider the case B < 0. The multiplicity of the origin is n — 2. There are 
two non-zero periodic solutions. In a neighborhood c = 0, q(c) = Bc n ~ 2 + 0(c™ -1 ). 
If n is odd then g(c) has the sign of c if |c| is large. Therefore, q has a positive 
solution and a negative solution if it is defined for large |c| . If n is even then q(c) > 
if |c| is large. Again, q has a positive solution and a negative solution if it is defined 
for large |c|. 

In the case B = and A ^ 0, the multiplicity of the origin is n — 1. There is 
only one non-zero solution which is a real solution because complex solutions occur 
in conjugates pairs. The argument used in the first case implies that this solution 
is positive if A < and is negative if A > 0. 

If we start with A = B = 0, then the multiplicity of z = is n. Perturb the 
equation, so that B = and A ^ 0. Since the total number of periodic solutions is 
unchanged by small perturbations, a real solution will bifurcate out of the origin. 
Now, we make a second perturbation so that B is negative and such that \B\ is 
small compared to |^4|. A second real nonzero solution will bifurcate out of the 
origin. If the equation is perturbed such that A — but B is negative and small, 
then the stability will be reversed and two non-zero real periodic solutions bifurcate 
out of the origin; one is positive and one is negative. □ 
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4. RIGID SYSTEMS 

Consider the system 

x = Xx - y + x(R n -x(x,y) + R n - 2 (x,y) -\ \- Ri(x,y)) 

(4.1) 

y = x + Xy + y(R n ^i(x,y) + R n - 2 {x,y) H h Ri{x,y)), 

where Ri is a homogeneous polynomial of degree i. The system in polar coordinates 
becomes 

f = ^Rn-iicosO^iaB) + r n ~ 1 R n ^ 2 {cos9,sin9) H h r 2 i?i(cos6»,sin6>) + Ar 

9 = 1. 

This system is called a rigid system because the derivative of the angular variable is 
constant. It was shown in that if the origin is a center then J^* Rk (cos 9, sin 9)d9 = 
for all 1 < k < n — 1. In fact, these definite integrals are the first focal values of 
the system. It is clear that the origin is the only critical point and if it is a center 
then it is a uniformly isochronous center; see [Hj- Limit cycles of Ij4.1|l correspond 
to positive 2it— periodic solutions of 

J- = R n ^r n + Rn^r"- 1 + ■■■ + R ir 2 + Ar 
d9 

In |7| , the family of rigid systems 

x = Xx - y + xRn-i(x, y), 

y = x + Xy + yRn-i{x, y) 

was considered. It was shown that if A J Q 27r i?„_i(cos 9, sin 9)d9 < 0, then there is 
at most one limit cycle. Cubic rigid systems were considered in 8 . It was shown 
that there are at most two limit cycles. 
Now, we consider the rigid system 

x = —y + x(Rn-i(x,y) + R„,-2(x,y) + Rn-3(x,y)) 

(4.2) 

y = x + y(Rn-i(x,y) + R n -2(x,y) + R n -z{x,y)). 
In polar coordinates, the system becomes 

^ l > ^Ji i d ~Ji— 1 , i ) — 2 

— = Rn-ir + Rn-2T + Rn-ZT 

d9 

If a function Ri does not change sign, then it is necessary to assume that i is 
even. We assume that n is odd. Hence, R n -2 is a homogeneous polynomial in 
sin# and cos# of odd degree; Therefore, i?„_2(cos 9, sin 9)d9 = 0. On the other 
hand, real periodic solutions occur in pairs. If ip(9) is a real periodic solution then 
so is — <p(9 + 7r). Let B = J Q 27r i?„_3(cos 9, sin 9)d9. The Liapunov quantities of 
the system are B and 1. If B = then the origin is unstable. By perturbing the 
coefficients of R n -3 such that B < and is small enough, the origin will be unstable 
and a small-amplitude limit cycle appears through multiple Hopf bifurcation. The 
following result follows directly from these remarks and Theorem 13.21 



Theorem 4.1. Suppose that n is odd, R n -i — 1, and R n -3 < 0. Let B 

1(1 



Jq V Rn-3 (cos 0, sin 9)d9. 



(i) The origin is not a center for system (|4.2[1 . 

(ii) If B < then system B4.2JI has at most one limit cycle. If B is small enough 
then the system has a unique limit cycle; this limit cycle is unstable. 
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(iii) If B = then system Ij4.2|l does not have a limit cycle. 

Remark 4.2. If the leading coefficient Rn-i does not vanish anywhere then the 
transformation of the independent variable 



8 i— > exp( / i? n _i(cosu, smu)du) 
J o 

reduces the polar equation into a similar equation but with a leading coefficient 
equals one. 
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